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Disclosed Invention: 

Device and Method for Improved Imaging in Nuclear Medicine 

Nuclear medicine imaging represents a challenge since the energies utilized 
are relatively high in comparison to Diagnostic Radiology requirements 
and the radiation source (radionuclides) is not directional and its 
distrition within the body is not precisely known. An extreme example is 
PET (Positron Emission Tomography) which involves very high energies. 
Emitted radiation which wscapes the body can be scattered which can 
alter its direction and energy. It is desirable (in most cases) to reject 
scatter radiation reaching the detector based on energy and/or direction. 

Most imaging systems used in Nuclear Medicine imaging are based on the 
Gamma camera design (a large scintillation crystal or array of scintillation 
crystals optically coupled to an aray of Photo Multuplier Tubes (PMTs) . 
Depending upon the imaging task, a flat or focused collimator is located 
between the Gamma camera and the object. This system provides a limited 
degree of spatial resolution and energy resolution while removing some 
fraction of the scatter radiation which would otherwise degrade image 
quality and diagnosis. Current systems use 1, 2, or 3 Gamma cameras detector 
units. At least one commercial system eliminates the use of scintillator 
crystals and PMTs with a semiconductor detector (CdZnTe) manufactured 
into a 2-D array. Drawbacks include the difficulty of growing large area, 
thick CdZnTe crystals which have problems with defects and forming a 
low noise, 2-D readout structure on a CdZnTe crystal. Grid collimation 
is still desirable for many applications .( 'DigiRad' , San Diego, CA 92126, 
858-578-5300) . Yet another design which avoids the use of a collimator 
(basically replacing it with a thin semiconductor which plays the role 
of a Compton scatter) can be thought of as a Compton Gamma camera. This 
design is not used commercially. 

A new Gamma camera design utilizes the concept of edge-on detectors as 
described by Nelson in patent 4,937,453. Manufacturing requirements in 
comparison to large area, thick, 2-D semiconductor arrays from CdZnTe, CdTe, 
GaAs, Ge, Si, SiC, HgI2, etc. are greatly reduced. Because radiation is 
incident on the edge of the detector the thickness of the detector now 
represents the maximum entrance aperture while the length or width of 
the detector area now represents the maximum attentuation distance for 
radiation. Drift and single-sided or double-side strip detectors 
configurations can be implemented depending on the semiconductor material 
chosen. Strip widths can be tapered (or curved in the. case of drift 
detectors) if focusing is desired. In the case of double-sided strip 
detectors electrons and holes can be collected to provide 2-D position 
information across the aperture or depth of interaction information can 



be acquired. 



Nelson has shown in patent 4,937,453 how large area 2-D detecors can be 
formed by the close proximity of multiple edge-on linear detector arrays. 
(See FIG. 1.) Limited focusing, in the case of strip detectors, is 
possible by tapering the strips for each linear array. Of course, a slat 
or focused grid can still be used with this design. Alternative geometric 
configurations may also be implemented other than a planar array (a 
conventional format) . FIG. 2a and FIG. 3 show arrays of edge-on detectors 
used in focused cylindrical and spherical lens-like geometries. It is 
obvious that these configurations can be "f illed-out " to form complete 
cylinders (rings) or nearly-complete spherical detector arrays. Collimation 
can be provided for individual edge-on detectors (FIG. 2b) and sheilding can 
be used as needed to minimize radiation escaping from neighboring detectors. 
The entire detector array mount can be moved as can individual edge -on 
detectors within the array. This allows the detector array to be refocused 
as needed or to reposition individual edge-on detectors in order to optimize 
its field of view for the specific distribution of radionuclide. The entire 
array or individual edge -on detectors can be dithered to provide appropriate 
sampling of spatial regions which would otherwise be 'dead areas' due to a 
lack of detectors at those postions. Alternatives to conventional collimation 
(which may also be used to expand the aperture of an edge -on detector and 
energy resolution) are the use of one or more of reflective x-ray mirrors 
or Bragg crystals, capillary x-ray optics, or refractive x-ray lenses. A 
combination of a x-ray mirror or Bragg crystal with a minifying capillary 
x-ray optic is shown in FIG. 4. The use of such devices must take into 
consideration the energy of the radiation and the physical size limitations 
of the (edge -on) Gamma Camera. 

Although this device and method have been described in terms of use for 
Nuclear Medicine, it is clear that industrial/scientific applications could 
also benefit. In addition, radiological applications such as Compton scatter 
imaging will also benefit from this design. 

Cooling of the edge -on detectors may be implemented as needed. Preferably 
the detectors will operate at or near room temperature for applications 
in Nuclear Medicine, but this may not always be possible. In some instances 
heating may be needed depending upon the detector package and the operational 
environement . 



October 17, 1999 



Robert S. Nelson 
2922 Upshur Street 
San Diego, CA 92106 



Disclosed Invention : 

Device and Method for Improved Scanned Mammography 



Developement goals for x-ray mammogrpahy imaging systems include high 
detectivity of disease conditions with minimal acceptable risk to the 
patient. Numerous efforts have been made or are underway to replace standard 
x-ray film-screen detectors with digital detectors. These digital detectors 
can have geometries which vary from slits and slots (which provide a level 
of self -collimation) to flate plates or 2-D detector arrays (which 
typically use a grid depending on the size of the breast being imaged) . 
Slit and slot devices are scanned across a stationary, compressed breast 
and thereby acquiring a single continuous image (the same format which 
is acquired when using film-screen or plate/2-D digitial detectors) . 

Advantages which digital detectors may hold over analog fim- screen 
detectors include an extended dynamic range, superior detection efficiency, 
and the option of easy use of post-processing of the image. In some 
instances the contrast resolution of the digital detector may exceed that 
of the film-screen unit over the frequency range of interest. 

Other aspects of the imaging chain can be optimized to further improve the 
final image. For example, the xray beam spectrum in a slit or slot format 
can be filtered using a variety of techniques (capillary optics, refractive 
lenses, multilayer mirrors or Bragg crystal spectrometers) or combinations 
thereof (see FIG. 3) . If the detector itself is capable of providing 
energy resolution, this could be used to improve contrast or reject a 
fraction of the incident scattered x-rays. The variable magnification from 
center-to-edge of the beam could be corrected with appropriate x-ray 
optics. This would reduce the problem of depth-dependent magnf ication for 
objects positioned at different locations within the breast. In addition, 
the x-ray optics may be used to focus the x-ray beam, thus increasing the 
intensity for a slit or slot scan format and thereby reducing the power 
requirement for the x-ray tube. 

Slit or slot scanning techniques are typically much more efficient at 
reducing the level of scattered x-rays reaching the detector (which 
reduces contrast resolution) . If slit or slot scanning is implemented such 
that the scanning speed minimizes the effects of patient motion (blurring) 
on contrast resolution of the local image then there may not be a need to 
acquire the entire image in a single, continuous pass of the scanning 
system. This suggests that greater local compression of the breast can be 
acceptable. The benefits include reduced x-ray tube heat loading (lower 
instantaneous power requirements) , reduced scatter, reduced x-ray attenuatio 
and filtration. Greater local compression can be attained using compression 



plates which are reduced in size. This can be attained using flat, 
contoured, or sliding compression plates (See Nelson and Zach, patent serial 
number 08/597,447 and pending CIP) as shown in FIG. la, lb. Image 
acquisition then involves acquiring a number of sub-images (with suitable 
overlap) which can be viewed separately or assembled into a larger image 
(see FIG. 2) . Overall scanning times will increase, but this should not be 
an issue since image acquisition times for current systems are very short. 
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Title of Invention: Device and Method for Improved Compton Imaging 

Date: June 2, 2002 

Description: 

Compton cameras typically consist of a Compton scattering device such as a thin 
semiconductor array aligned with (and typically separate from) a thicker detector. 
A pending application tilted "Device and System for Improved Imaging in Nuclear 
Medicine and Mammography" describes an implementation that includes such a thin 
array used with an edge-on or strip detector array -based Gamma Camera. This is a 
traditional format for a Compton Gamma camera. 

Now consider that our edge-on (or strip) detector has sufficient spatial resolution so that 
an event (photoabsorption or scattering) occurring within a detector element can be 
localized. This can be achieved by making the edge-on detector into a pixilated array or 
by having crossed-readout stripes on opposite sides of the edge-on detector. Migration 
times for holes and electrons can be used to localize the event with respect to the two 
sides. There are a number of configurations such that the edge-on or strip detector array 
can function as a Compton camera: 

1. The semiconductor edge-on (or strip) detectors act as both Compton scatterers 
and Compton scatter detectors. (See Fig. 1.) 

2. The semiconductor edge-on (or strip) detectors are constructed from a stack of 
edge-on (or strip) detectors made from different materials. One material may be 
selected as the Compton scatter based on a high likelihood that for appropriate 
photon energies scattering dominates absorption. The second material may be 
selected such that absorption dominates scatter. We are not limited to using only 
two materials, the idea is easily extended to use 3,4,... materials. (See Fig. 2.) 

3. Alternate semiconductor edge-on (or strip) detector materials within an array 
configuration. The materials can be selected using the criteria described in 
configuration 2. (See Fig. 3.) 

Configurations 2 and 3 may be more complex to implement than configuration 1 but may 
improve Compton imaging efficiency relative to configuration 1. All 3 configurations 
may be used for conventional imaging. 
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Model calculations of the response of CZT strip detectors 
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ABSTRACT 

Position-sensitive CZT detectors for research in astrophysics in the five - several hundred keV range are being 
developed by several groups. These are very promising for large area detector arrays in coded mask imagers and 
small-area focal plane detectors for focusing X-ray telescopes. We have developed detectors with crossed-strip readout 
and optimized strip widths and gaps to improve energy resolution. A "steering electrode" is employed between the 
anode strips to improve charge collection. A model, of charge drift in the detectors and charge induction on the 
electrodes has been developed to allow us to better understand these types of detectors and improve their design. 
The model presently accounts for the electric field within the detector, the charges' trajectories, mobility and trapping 
of holes and electrons, and charge induction on all electrodes including their time dependence. Additional effects 
are being added. The model is described and its predictions are compared with laboratory measurements. Results 
include (1) the dependence of anode, cathode and steering electrode signals on interaction depth, transverse (to the 
strips) position, electron and hole trapping, strip width and gap, and bias, (2) trajectories of charges for various 
anode and steering electrode bias voltages, (3) a method to improve energy resolution by using depth of interaction 
information, and (4) an electrode geometry and bias optimized for the improved energy resolution. In general, the 
model provides good agreement with the measurements. 

Keywords: CdZnTe, CZT detectors, strip detectors, solid state detectors, X-ray astronomy 

1. INTRODUCTION 

CdZnTe (CZT), a semiconductor detector of X-rays and gamma-rays, has numerous applications in scientific and 
technological fields, most notably in medical imaging and high energy astrophysics. An important objective of high 
energy X-ray astronomy is a full sky survey at energies — 10 to several hundred keV at a sensitivity level that will 
yield many new sources. 1 These sources include ordinary objects such as stellar coronea, galaxies, cluster of galaxies 
as well as exotic ones, such as black holes, active galactic nuclei and neutron stars. Fundamental information on many 
topics in astrophysics is obtained by observing those sources in X-rays. Coded mask techniques can be used along 
with a detector having submillimeter position resolution and a few keV energy resolution to achieve this objective. 2 
(For more information on the objectives and status of high energy astronomy, see Matteson et al. 3 ) 

These requirements are met with CZT. It is dominated by photoelectric absorption due to high Z, and its large 
bandgap (~ 1.4 eV) allows room temperature operation. Its high bulk resistivity at room temperatures, ~ 10 11 fkm, 
results in low leakage currents and noise in a small detector, i.e. 1 cm 2 . CZT's energy response is linear over the energy 
range of interest. 4 However, CZT material has poor hole transport properties. Conventional planar configuration 
for the detector shows low energy tailing in their spectra for energies higher than 60 keV due to incomplete hole 
collection. Strip or pixel detectors can eliminate this problem by virtue of the small pixel effect. 5,6 

Further author information: 
E. Kalemci: E-mail: emrahk@mamacass.ucsd.edu 
J. L. Matteson: E-mail: jmatteson@ucsd.edu 
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Figure 1. Line profiles and energy resolution of prototype CZT strip detector 

Position-sensitive strip CZT X-Ray detectors for use in X-Ray Astronomy have been developed by UCSD and 
WU which have very good energy resolution (e.g. 4.2 keV at 60 keV and 6.4 keV at 122 keV). (See Fig. 1) With 
lower noise electronics, such as described in Matteson et al. 7 , even better resolution can be obtained. 

In order to obtain a better understanding of the performance of the detector, a modelling program was developed. 
■f.^j. , The objective is to predict the spectra and the signal waveforms of X-rays interacting in a CZT strip detector. This 
fgjm* requires determining the positions of interactions, propagating the electron-hole pairs, calculating the charges induced 
™- • on electrodes and modelling pulse shaping electronics. At this time, the model does not account for X-ray interactions 

* f and electronic effects, but its calculations of charge transport and signal production are in very good agreement with 

* ■ the experimental measurements. 

2. BASIC PRINCIPLES FOR THE MODEL 

2.1. X-Ray Interactions 

For the energy range of interest, ~ five - several hundred keV, interactions of photons with the detector are mostly 
photoelectric. For example at 100 keV photoelectric cross section is ~ 10 times greater than the Compton scattering 
cross section. At 250 keV, the cross sections are equal. 

Cd and Te have K edges at ~ 25 keV and ~ 31 keV respectively. So an X-ray with energy greater than 25 keV 
can produce a photoelectron and a characteristic K X-ray. Normally, the K X-ray will be absorbed within ~ 100 fim 
of the interaction position, unless the interaction is close to the physical boundaries of the detector, and it escapes. 
The photoelectron loses energy by creating electron-hole pairs. It will be stopped completely within a short distance, 
e.g., 30 //m for 100 keV X-rays This, along with thermal diffusion, causes an initial spreading of the electron-hole 
plasma. 

2.2. Charge Transport 

When bias is applied between anodes and cathodes, charges drift along the electric field lines. Up to fields of 
~ lOkV/cm the drift velocity of electrons and holes are proportional to the field. 8 (v — \i E> where fj, is the 
mobility.) For higher fields, they drift with a constant speed. In CZT, electrons move much faster than holes and 
typical mobilities are fi e ~ 1000 cm 2 V^s" 1 , ^ ~ 50 cm 2 V 1 s~ 1 , respectively. 9,10 A lower value of hole mobility, 
fj, h ~ 15cm 2 .V"" 1 s~ 1 , is also reported. 11 
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2-3. Trapping and Det rapping 

The charges can be trapped on their paths due to trapping sites in the crystal structure. The trapping can be 
quantified by trapping time r defined as mean free drift time of carriers before trapping occurs. 12 Both electrons 
and holes detrap as well, which restores detector's charge neutrality. A useful quantity is the trapping length defined 
as: L(x) = (/i r) E(x). If the electric field is constant, the charge remaining as a function of distance is written as: 

Qe,h(x) = 3^(0) exp[-d/LeAx)l (1) 

where Q(0) is the total number of electron-hole pairs created and d is the distance from the interaction site. Since 
holes move much slower than electrons, they are more susceptible to trapping. The typical values found in the 
literature for electron and hole /zr products are fx e r e ~ 3 x 10" 3 cm 2 /V, ii h r h ~ 1 x 1(T 4 to 8 x 10" 6 cm 2 / V. 10 » 13 
Those values correspond to a r e of around 3 /is and a r h of 2 /*s to 160 ns respectively. For a typical field of 1000 V /cm, 
the corresponding trapping lengths are L e ~ 3 cm and L h ~ 1 mm to 0.08 mm respectively. 

Detrapping is a complicated phenomenon whose effects depend on temperature, material and pulse shaping 
electronics. Its details are discussed by Martini et al. 14 and by Mayer. 15 

2.4. Charge Induction 

The charge induced on each electrode can be calculated by the weighting potential method, 16 Q ind = -Q(x) W(x), 
where W{x), the weighting potential, is a number between 0 and 1, and can be obtained by solving the Laplace 
equation in the system where unit potential is assigned to the particular electrode and the remaining electrodes are 
assumed to be grounded. Electron-hole pairs separate under the applied field, and, as they drift, they induce charges 
on each electrode with the relation above. Because of trapping, some charges do not reach the electrodes, and the 
total induced charge on an electrode is a function of an interaction's position. For a planar detector, an analytic 
solution to total charge induction as a function of depth is given by the Hecht equation 17 : 

Qind/Q(0) = (L e /D)(l - exp[-(D - x)/L e )) + (L h /D)(l - exp[-x/L h ]) (2) 

where D is the detector thickness and the interaction occurs at a distance x from the cathode. Note that the 
Hecht equation is valid only when detrapping within the integration time is negligible. For a strip detector the 
induced charge can be found numerically as in this paper, or can be calculated analytically by integrating the current 
equation. 5 * 18 

For a strip or pixel detector, an anode has strong weighting potential only in its vicinity due to its small size 
compared to the detector size. See Fig. 7 for an example. If illumination is from the cathode side, the holes trapped 
near the cathodes do not affect the anode signal since there is no appreciable weighting potential. Indeed, for 
the energy range in interest, most of the interactions are near the cathodes. Since electrons show good transport 
properties, they are collected at anodes. This is "the small pixel effect". Due to this effect, strip or pixel detectors 
yield superior performance to planar detectors. 5 ' 6 

3. MODELLING 

By using the principles explained above, . a computer simulation of a CZT detector was developed that predicts in- 
duced charge on each electrode for various geometries. It consists of two parts: (1) calculation of electric fields and 
weighting potentials, and (2) charge drift, trapping and calculations of induced charge using the fields and potentials. 
The former is done with a commercial program, MAXWELL, 19 which, for a given geometry calculates the potentials 
and fields and saves the values to a file. Then a C program reads those values, and drifts charges on a rectangular 
grid. The grid is coarse at the cathode side (10 /mi) and fine (1 /mi) near the anode side to accurately follow the 
trajectory. Cathodes are at the top and the electrons move down on grid points following the electric field lines. If 
the grid size is M, depth on the grid is z and the interaction depth is z Q , then for electrons: 



• Az = M 
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• Arr = /i c E x At 



• L e =f, e T e (EZ + E 2 2 ) 1 / 2 

• Q(z) = Q(z - 1) exp(-(x 2 4- M 2 fl 2 /L e ) and Q T (z - 1) = - Q(z - 1) 

• Qind = Q(z) W(z) + Q T (i) (i) 

where Q iTl d is the induced charge on a specific electrode and Q T {%) is trapped charge at site i. A similar set of 
equations applies to holes. Induced charge is calculated for each electrode type. 

A cautionary remark is needed here since for some geometries involving steering electrodes, E z might change sign 
close to steering electrodes. Those cases must be treated carefully. (For example, for an interaction in this region, 
the holes prefer to drift to steering electrodes, instead of cathodes.) 

The model assumes a linear field-velocity relation throughout the detector. The charges trapped are not de- 
trapped. Point charges are used, and there is no diffusion or initial spreading of electron hole plasma. The model 
does' not include Compton scattering and K X-rays. 



4. EXPERIMENTAL SET UP 

The detector modelled is a laboratory prototype which was developed to study techniques for X-ray imaging instru- 
ments such as HEXIS 7 and MARGIE. 20 There are 22 anode and 22 cathode orthogonal strips on the 12 x 12 x 2 mm 3 
detector. (See Fig. 2.) In addition, there are steering electrodes in between anodes to enhance the charge collection. 
Steering electrodes are connected to each other. Anodes and steering electrodes are 100 yum wide, cathodes are 
450 fxm wide, and all have 500 jzm pitch. The anodes are biased at 200 V and steering electrodes at 180 V, so 
charges created above steering electrodes are directed to anodes. (See Fig. 3) A ceramic carrier holds the detector. 
The electrodes are coupled to Amptek A250 charge sensitive preamplifiers. Signals then are processed by shaping, 
amplifying and triggering circuits, and digitized by ADC's which give the total signal. Ah on-line program processes 
the data and builds multiple spectra according to various event selection criteria. 7,21 Off-line processing is also used. 

The system also carries a planar detector with the same size as the strip detector, monitored by the same 
electronics. This gives an opportunity to compare the performance of strip detector to the planar detector. 

An additional setup was used to measure pulse shapes from various electrodes. The detector was illuminated by 
X-ray sources with two Amptek A250 preamplifiers connected. Signals were recorded using a Tektronix DS A 601 four 
channel digital oscilloscope. The archive of those measurements includes various combinations of anodes, cathodes 
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Figure 2. Prototype CZT strip detector. 22 anode 
strips with 500 microns pitch and 22 steering elec- 
trodes. Total size is 12x12x2 mm 3 . 



Figure 3. Detector electrode and bias network con- 
nections. 
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and steering electrodes. Sources used were 57 Co which primarily emits 122 keV 7-rays which are fully absorbed, and 
137 Cs which emits 7-rays at 662 keV. These can Compton scatter and leave energies up to 382 keV in the detector. 

To obtain transport properties, a set of measurements was made with an alpha particle source on a planar 
detector. The alpha beam was collimated to make it monoenergetic. The anode was coupled to an Amptek A250 
preamplifier and the detector was illuminated from cathode side to determine electron transport properties. Then 
the bias was reversed to get hole transport properties. For bias voltages from 15 V to 210 V, the signals were recorded 
with the digital oscilloscope and the spectra were recorded by a ND76 Spectrum Analyzer. 

5. RESULTS 

5.1. TRANSPORT PROPERTIES 

The measurements done using the planar detector gave important insights. 
Spectra obtained from the. motion of electrons had good energy resolution 
(~ 3% FWHM) for the electric fields greater than 1000 V/cm. The am- 
plitudes were used to fit the Hecht equation (See Fig. 4), and a (fir) e of 
5.5xl0~ 6 cm s V -1 was obtained. For the holes however, 3 (as shaping time 
of the amplifier caused a severe ballistic deficit effect. The peaks were very 
broad, ~80% at any bias. This may be due to non-uniform hole trapping 
and electric field. 8 

The individual signal versus time curves were also used to obtain trans- 
port parameters. The mobility of electrons was measured by using the 
slope of the linear part of signal vs time. 14,15 See Fig. 5 for an exam- 
ple. The values were between 900 and 1100 cm 2 V -1 s -1 with a mean 
of 960 cm 2 V -1 s _1 . But even for high biases, signal versus time curves 
were not completely linear, and showed a gradual non-linear increase to the 
maximum. This is a typical sign of detrapping 14,15 which makes the Hecht 
equation fit questionable. Still, it gives an effective jut, and determined 
value obtained, 5.5 x 10~ 6 cmsV" 1 , is close to 6.2 x 10~ 6 cms V" 1 which 
resulted in the best agreement with the measurements. (Fig. 13, Sect. 5.3.) 

The slope method was also applied for the holes, and an average of 
fXh ~ 22 cm 2 V^ 1 s -1 was obtained, which is smaller than usually found in the literature. Perhaps, an effective 
mobility was found with this method, since hole transit time (time required for carriers to tran verse the detector) 
may be comparable to trapping and detrapping times. 15 

5.2. TRAJECTORIES AND WEIGHTING POTENTIALS 
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Figure 4. Measured peak height as a 
function of applied bias for alpha par- 
ticles impinging on the cathode. Solid 
line is a fit to Hecht Equation, eqn(2). 




Figure 5. Signal versus time of an electrode when an alpha particle impinges on the cathode for a planar detector. 
The x-axis unit is 100 ns per division, and y- axis unit is 20 mV per division. 
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Figure 6. Model calculation of charge drift trajectories. (X- and y-axes are not scaled the same.) Filled boxes are 
anodes and empty boxes are steering electrodes here and in Fig. 7 
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Figure 7. (a) Model calculation of weighting potential of the prototype detector. The dashed lines are steering 
electrodes weighting potential, (b) Model calculation of weighting potential without the steering electrode 

Results of electric field calculations are shown in Fig. 6. It is clearly seen that 
most of the trajectories are directed onto the anodes by the steering electrode. 
However, some go to the steering electrodes and the gaps. For interactions 
at random positons near to the cathodes, 85% go to the anodes, 10% to the 
steering electrode, and 5% to the gap. 

Figure 7a shows weighting potential distributions of all three types of elec- 
trodes. The cathode weighting potential extends deep into the detector, e.g., 
W ~ 0.1 at a depth of 1.2 mm. Steering electrode weighting potential is also 
very extended, W - 0.1 at depth 0.4 mm. As discussed before, if weighting 
potential of an electrode is extended, signals from that electrode are affected 
by hole trapping. On the other hand, the weighting potential of the anodes 
is localized to the immediate vicinity of the anodes. Thus, the anodes benefit 
from the small pixel effect: 

It was well established that addition of steering electrodes to the detector 
improved energy resolution. 2 Figure 7b explains the main reason for this. 
Without the steering electrodes, weighting potential of the anode extends fur- 
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Figure 8. Effects of the steering 
electrode and nearest neighbor sum 
for an anode spectrum. Source is 
241 Am. From Matteson et al. 2 
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ther, reducing the small pixel effect. The effect of relative biasing is also modelled. When there is no bias between 
steering electrodes and anodes, 20% of the trajectories end tip in the gap, which would result in incomplete charge 
collection. Also from the remaining 80%, half of them end up in steering electrode which does not provide position 
information. Fig. 8 shows the difference in spectra when there is biasing and there is no biasing. 

5.3. DEPTH OP INTERACTION 

Weighting potential curves clearly indicate that the amplitude of the signals is a function of depth of interaction. This 
effect, which causes a low energy tail due to incomplete charge collection, was studied with the model by calculating 
the total signals for interactions at various depths. An electron trapping length of 6 cm and hole trapping length 
of 0.24 mm were used. Fig. 9 shows the results. For the anodes, the response is uniform to ±2% for depths less 
than 1.7 mm. For the depth greater than 1.8 mm, the anode signal amplitude becomes lateral position dependent. 
Typically, useful signals are produced on several neighboring cathode strips due to trapping of holes. Fig. 10 is an 
example of this which shows signals of two neighboring cathodes as a function of time. For the single cathode with 
the largest signal, Fig. 9 shows a rapid decrease in signal with depth. When the two neighboring cathode signals are 
added, the depth dependence is much less severe. Therefore by adding neighboring cathode signals to the one which 
has the maximum signal, better spectra can be obtained (See Fig.ll). 

The model was tested with a series of experiments. Data 
were collected by flood illumination with 57 Co at 122 keV. For 
each event, the ratio of cathode to anode signal and steering 
electrode signal were determined. These were plotted as a 
scatter plot in Figure 12 along with the model predicitons. 
The agreement is very good. This technique was used by Slavis 
et al. 21 to reduce background by rejecting events which were 
not incident X-rays. 

It has been mentioned before that cathode to anode ra- 
tio is a monotonic function of depth. 2 Relative anode signal 
versus cathode to anode ratio is plotted along with the model 
prediction. The agreement is very good, except at a very small 
region with high cathode to anode ratios. (See Figure 13.) 

The model results were used to correct for the dependence 
of signal on depth of interaction and thus improve energy res- 
olution. Pulse heights were multiplied by correction factors 
which were obtained from the model, using the cathode to an- 
ode ratio. Results are shown in Figure 14. In the corrected 
spectrum, interactions deeper than 1.7 mm have been dis- 
carded. This means discarding 18% of the events for 100 keV 
and less than 1% for 60 keV. There is significant improvement 
in energy resolution. The corrected spectrum has energy resolution of 7 keV whereas the uncorrected neighbor 
summed spectrum has 10.5 keV energy resolution at 122 keV. For comparison, a spectrum obtained from the top 
0.05 mm of the detector is also plotted. This region presumably has the best energy resolution. 22,2 There is almost 
no difference between two spectra. The low energy tail is eliminated and the shape is symmetric. This method can 
be implemented in the software very easily. 
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Figure 9. Model calculations of total signals at 
electrodes of prototype strip detector and planar de- 
tector. 



5.4. SIGNAL VERSUS TIME 

Depending on the interaction position, there are various types of signal behavior of the electrodes. Several examples 
are shown in the following, which are modelled with following parameters: r e = 6 /xs , \x t — 1000 cm 2 V*" 1 s -1 and 
r h ~ 645ns and ^ = 40cm 2 V" 1 s -1 . Fig.15 shows the general results for a typical interaction at 200 /zm. The anode 
signal rises slowly until electrons reach significant weighting potential (See Fig. 7a). Then, it rises rapidly with a 
typical risetime of 50-70 ns. Only 3% of the electrons are trapped. On the other hand, the cathode signal increases 
promptly due to its extended weighting potential. This signal usually consists of two parts, (1) a rapid component 
due to electrons drifting away from the cathode (they move to lower weighting potential) and (2) a slow component 
due to holes approaching the cathodes. Charge collection is incomplete due to trapping. Fig. 15e shows a typical 
signal with all the features described above, measured with our setup with 57 Co. 
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Figure 10. Typical cathode and neighboring cathode 
signals as a function of time. X-axis is 100 ns per 
division and y-axis is 1 mV per division. 



Figure 11. Effect of nearest neighbor sum for a cath- 
ode spectrum. Source is Am 241 . From Matteson et 
al. 2 
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Figure 12. Modelled (left) and measured (right) 
depth of interaction effect on cathode to anode ratio 
and steering electrode signal for 122 keV X-rays. 




Figure 13. Anode signal versus cathode to anode 
ratio. Triangles are measurements for 122 keV and 
solid line is the model prediction. 



A second example (Figure 16) illustrates the behavior of signals on the steering electrode when electrons are 
collected at an anode. The interaction is at 600 pm. The signal first rises rapidly but then drops to a negative value. 
Then it rises very slowly. To understand this, recall Fig. 7a. The weighting potential of the steering electrodes rises, 
down to a depth of 1.6 mm, then drops quickly to zero near the anodes. The modelled signal drops down not to 
zero but negative values because of trapped holes. The slow increase at the end is due to untrapped holes going 
away from the steering electrodes. The measurements (See Fig. 16e for a typical signal) show those properties. For 
the case of an event at the very top, like the first example, steering electrode signal drops to zero (Figure 15b) since 
most of the holes are not trapped and weighting potential is very low. One discrepancy between simulations and 
measurements is the fall time of the steering electrode signal after electrons are collected by anodes. Fall times of 
~ 50 ns are predicted, but instead, ~ 150 ns are observed. 

Results of electric field calculations show that for some of the interactions, the electrons can drift to the steering 
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Figure 14. Comparison of depth corrected spectrum (histogram) with uncorrected spectra. Neighbor summed 
single spectrum is obtained from the top 0.05 mm of the detector. See Section 5.3 for details. 

electrodes, or even to the gaps between steering electrodes and anodes. These two cases were modelled also. Fig. 
17 shows the case where electrons end up at a steering electrode for an interaction at 0.2 mm. The rise time on the 
steering electrode signal is greater than that on the anode. The anode signal in this case, rises first and goes down 
to almost zero since weighting potential is very low at the depths where trapped holes are located. Fig. 17e shows a 
measured steering electrode signal for Cs 137 which agrees with the model. 

What happens when electrons reach the gaps is unknown. It is modelled with idealized surface conditions for 
which an electron stops when it reaches the surface. The surface of a semiconductor crystal has a lot of imperfections 
and therefore has a lot-of trap centers for the free charges. Basically, the model treats electrons as. trapped when 
they reach the surface. Fig. 18 shows an example. Both anode and steering electrode do not collect full charge. The 
steering electrode signal (Fig. 18b) first increases to the maximum as usual, then decrease some, since electrons are 
directed to the anodes. However, they can not reach the anode, therefore induce a positive signal at the steering 
electrode. (Compare with the case when electrons reach the anode, the steering electrode signal is negative in this case 
due to trapped holes.) One of the measurements for the steering electrode showed a similar behavior expected from 
such interaction (Figure 18e). The timescales and amplitude are consistent with the model. If this interpretation 
is true, then the electrons move slower very close to the surface and then trapped. The slow increase at greater 
timescales is due to holes moving away from the steering electrodes. 

6. CONCLUSIONS AND FUTURE DIRECTIONS 

The model is able to predict correctly the fundamentals of how our detector works. In terms of predicting total 
signals, it agrees within 1% with the measurements (See Fig. 13). This shows that some parameters that are not 
included in the program, such as diffusion, and range of photoelectron do not affect the performance of the detector 
very much. In terms of temporal behavior, the model is within 10% -20% when parameters that are usually found in 
literature are used. The discrepancies may be due to pulse shaping of electronics or detrapping which was observed 
for the planar detector. (See section 5.1.) 

Trajectory and weighting potential calculations helped us to understand the improvement obtained with adding 
steering electrode to the system. It is also shown that, the small pixel effect is a result of not only the size of the 
electrode, but also the distance between electrodes (See Fig. 7). With this in mind, an optimized geometery for 
enhancing the small pixel effect was studied. In this geometry, anodes are 50 fxm wide and steering electrodes are 
250 /im wide. Pitch size is 500 /xm. More bias is applied between steering electrodes and anodes in order to decrease 
trajectories ending at wide steering electrodes. Preliminary model results showed an increase from 85% to > 90% in 
relative anode signal for very deep interactions compared to detector used now. A detector with this design will be 
tested later this year. 
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Figure 15. Modelled, (a)-(c), and measured waveforms, (e), and modelled trajectory, (d), for an interaction at 
200 /mi depth. In (e), the cathode (C) gain is 1.67 times the anode (A) gain. The x-axis unit is 100 ns per division 
and y-axis unit is 1 mV per division. The overshoot and ringing of anode signal is due to preamlifier response to the 
short risetime of the signal. See Sect. 5.4 for details. 

Both cathode and steering electrode signals are strongly depth dependent. Depth of interaction information, 
which can be obtained using the model, can be used to reject background 21 and to improve the energy resolution, 
as demonstrated in this paper. 

The model will be improved by adding a Monte-Carlo code to simulate X-ray interactions in the detector. 
Detrapping and modelling of electronics will also be included. Signals obtained from each interaction will be used to 
build a spectrum, so that the model will be able to predict the energy resolution of detectors with different geometries 
and biasing schemes for the electrodes. 

The predictions of trajectory calculations will be tested by localizing the X-ray interaction with a 30 collimated 
beam which is positioned by a precision X-Y stage. 4 In particular, this will allow us to study the case when electrons 
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Figure 16. Modelled, (a)-(c), and measured waveforms, (e), and modelled trajectory, (d), for an interaction at 
600 prn depth. In Fig. (e), the steering electrode (S) gain is 1.67 times the anode (A) gain. The x-axis unit is 100 
ns per division and y-axis unit is 2 mV per division. See Sect. 5.4 for details 



reach the gap. 

More experiments will be done on planar detector to understand the transport properties, especially detrapping. 
Homogeneity of the material will be studied by using the collimated beam. 
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Figure 17. Modelled, (a)-(c), and measured, (e), waveforms and trajectory, (d), for an interaction in which electrons 
drift to steering electrodes. The depth is 200 ^m. In (e), the x-axis is unit 100 ns per division and y-axis unit is 2 
mV per division. Events ending on the surface of steering electrodes constitue ~ 10% of all interactions. See Sect. 
5.4 for details. 
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Figure 18. Modelled waveforms, (a)-(c), and trajectory, (d), for an interaction at a depth of 700 /urn. In this case, 
electrons never reach the anode, but collect in the gap. Fig. 17e is a measured waveform which appears to be a gap 
event. The x-axis unit is 100 ns per division and y-axis unit is 2 mV per division. See Sect. 5.4 for an extended 
discussion. 
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Depth Measurement in a Germanium Strip Detector 

E. A. Wulf, J. Ampe, W. N. Johnson, R. A. Kroeger, J. D. Kurfess, and B. F. Phlips 



Abstract — We have demonstrated the ability to determine the 
depth of a gam ma- ray interaction point over the full active volume 
of a thick germanium strip detector This capability provides depth 
resolution of less than 0.5-mm full width half maximum (FWHM) 
at 122 keV in a device 11 mm thick with a 2 -mm strip pitch. Fifty 
channels of electronics have been developed and tested with a 25 x 
25 germanium orthogonal strip detectors. Experiments examining 
the capabilities of the system and demonstrating a simple Compton 
telescope using a single detector have been performed. 

Index Terms — Gamma-ray imaging, position measurement, 
radiation detectors, semiconductor devices, spectroscopy, time 
measurement 



I. Introduction 

GERMANIUM strip detectors combine excellent en- 
ergy resolution for gamma-ray detection with good 
two-dimensional (2-D) resolution. With the addition of depth 
information these detectors have excellent overall position 
resolution. Orthogonal strips on the front and rear faces of the 
crystal allow germanium strip detectors to locate a gamma-ray 
interaction in two dimensions accurate to the width of the 
strips. A gamma-ray interacts in the crystal and its position 
is determined by the intersection of the triggered strips on 
opposite sides of the detector [1]. The depth of the interaction 
is determined by looking at the timing difference between 
signals from collection of holes on one side of the detector 
and electrons on the other side as was recently demonstrated 
by [2] and [3]. The excellent energy resolution of germanium 
detectors makes it possible to determine if the signals collected 
at the front and back are from the same gamma-ray interaction. 
A germanium detector with submillimeter resolution in three 
dimensions is of interest in gamma-ray astrophysics for the next 
generation of instruments. It should also have the potential to 
improve the resolution of positron emission tomography (PET) 
[4]. Another application is for the GRETA detector under study 
for use in nuclear physics experiments by the Department of 
Energy [5]. 

The detector used for this work and the work of [2] is a 25 x 
25 germanium orthogonal strip detector with 0.2-cm strip pitch 
that is 5 x 5 x 1.1 cm deep [1]. It has lithium strips held at 
4-1.5-kV bias potential that collect electrons and boron strips 
on the opposite face to collect the holes. 
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The interaction depth is directly related to the time difference 
between when the electron and hole signals are collected on op- 
posite sides of the detector. This can be seen in Fig. 1 , which 
shows that for an event occurring near the boron face of the de- 
tector it takes 113 ns for the electrons to travel to the lithium 
face. The simplest way to conceptualize measuring the time dif- 
ference between charge collection is to determine the difference 
between the times when each of the preamplified signals crosses 
50% of its total value. The timing difference in a 1.1 -cm-thick 
detector is approximately ± 1 20 ns for conversion near the front 
or the back of the detector. 

II. Electronics 

To instrument a detector for depth information, one must de- 
termine the time difference between charge collection as well 
as the energy of the interaction. To determine the energy of an 
interaction, shaping amplifiers and analog to digital converters 
(ADC) are needed for all 50 strips. The depth determination 
requires a discriminator on each of the strips on the front and 
back of the detector and a time to digital converter (TDC) or the 
equivalent to measure the relative timing of the signals. 

A major design question is whether a constant fraction 
discriminator (CFD) is necessary or if a simple leading edge 
discriminator (LED) is adequate for the relative timing of the 
signal rise. One type of CFD works by making two copies of 
the input signal, inverting and delaying one copy, attenuating 
the amplitude of the other and adding the two signals. This 
creates a zero crossing that occurs when the original signal was 
a fixed percentage of its full value. This is useful for eliminating 
time walk as a function of amplitude. The problem with CFDs 
in this application is that 1 50 ns of delay are necessary. This 
may be difficult to implement in future compact, low -power 
electronics. In contrast, an LED triggers when the input signal 
goes over a specific voltage and therefore can trigger at different 
times for different pulse amplitudes. This may not be a large 
issue for the germanium strip detector because the amplitude 
of the signals on the front and back of the detector are the 
same and the time walk is expected to be similar. 

To read out all 50 strips on the detector with both energy and 
depth information requires 50 channels of shaping amplifiers, 
50 channels of discriminators, 50 channels of ADCs, and 50 
TDC channels. A decision was made to create a Nuclear Instru- 
ment Module (NIM) module that incorporated the shaping and 
discriminator functions in order to reduce the total number of 
modules. The outputs from eV5093 preamplifiers are fed into 
a buffer amplifier with a gain of 20 (see Fig. 2). The signal is 
then split and one copy is shaped by a four-pole shaper with a 
fixed gain and is fed to an ADC. Another copy of the amplified 
preamplifier signal goes to a fast shaper with an integration and 
differentiation time of 50 ns. The output of the fast shaper is run 
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Fig. 1. The digitized preamplified signals from a germanium strip detector. 
The dashed curve is the signal as holes are collected on the boron side of the 
detector which is closest to the 241 Am source. The solid curve is from the 
lithium side as the electrons are collected. There is a 1 1 3-ns difference in the 
time when the signals reach their midpoint [2]. 
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Fig. 2. A schematic of the NIM electronics constructed at NRL and used to 
instrument the 25 X 25 germanium strip detector. 




Fig. 3. A picture of one of the NRL electronics boards with four channels. 
Four of these boards are packaged together to produce one double-wide NIM 
module. The module to the left is a board holding multiple eV5093 preamplifiers 
that were used to produce test signals. 
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Fig. 4. Spectrum of 57 Co source as measured with the germanium strip 
detector and the NIM electronics. 



to a discriminator and compared to a dc level set by a front panel 
potentiometer. The output of the discriminator is summed with 
all other channels and used as a master trigger to start the ADC 
and as a common stop for the TDC. Another copy of the dis- 
criminator output is used to enable the comparator used in the 
CFD electronics. 

The CFD section is composed of two copies of the amplified 
preamplifier signal. One copy is delayed by 1 50 ns and the other 
is attenuated to 50% of its original amplitude. These two signals 
are fed into a comparator that fires when the two signals have 
the same amplitude. In effect, this produces a signal when the 
preamplifier signal has risen to half of its total value. The CFD 
signal is used to start a TDC channel for each front and back 
strip. 

Each electronics board supports four detector channels and 
four boards are included in one double wide NIM module (see 
Fig. 3). The outputs from these modules are fed into TDCs and 
ADCs residing in a CAMAC crate which is read out by a PC 
running Linux. The data is recorded on an event by event basis 
and saved to disk and tape for later analysis. This system main- 
tains the excellent energy resolution, 1.6 keV at 122 keV, of a 
germanium detector as can be seen from a typical 57 Co spec- 
trum in Fig. 4. 



III. Depth Measurements 
A. Detector Attenuation 

The depth capabilities of the detector are demonstrated by ob- 
serving the attenuation of gamma-rays as they pass through the 
detector. These tests confirmed that the depth of the interaction 
could be measured but are not an accurate way to determine 
the actual depth resolution of the system. The attenuation ex- 
periment was done by placing a source near the boron face of 
the detector and producing a histogram of the time difference 
in charge collection between the boron and lithium face. Each 
event histogramed had to have only one strip with a signal on 
both the lithium and boron side and each signal had to be the 
correct energy to within 5 keV. 

o7 Co has a 1 22 keV gamma-ray line that is attenuated 85% by 
the detector volume. The radiation length is 5.75 mm, which is 
approximately half the detector thickness. A plot of the number 
of counts as a function of time difference between charge col- 
lection on the front and back face is shown in Fig. 5. The face 
of the detector that was closest to the source was the boron face 
which corresponds to negative time differences and the lithium 
face to positive differences. The theoretical exponential atten- 
uation of the germanium is shown superimposed as the dashed 
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Fig. 5 . The number of photo-peak events for the 1 22-ke V gamma- ray line from 
5 7 Co as a function of time difference between charge collection on the boron and 
lithium face. The source was placed 40 cm from die boron face of the detector. 
The dashed line is die theoretical exponential attenuation of the gamma-rays by 
the germanium that makes up the detector. 



line on the plot. The total time difference is shown to be 215 ns 
for the 1.1 -cm-thick detector. This is similar to the total time 
difference found by [2]. The 15 ns difference between hole col- 
lection on the boron side and electron collection on the lithium 
side is due to the difference in drift velocities in germanium. At 
the detector's bias voltage of 1 500 V and a temperature of 80 K, 
the drift velocity of the holes is 7.5 x 10° cm/s, and that of the 
electrons is 83 x 10 6 cni/s [6]. Using a detector thickness of 
1.1 cm yields a total time difference that is larger than observed 
by 25%. This suggests a small nonlinearity in depth timing near 
the surfaces. The experiment was also performed with 241 Am 
and 137 Cs, which showed good agreement with the theoretical 
exponential attenuation curves. The differences between the at- 
tenuation curve and the measured values are most probably due 
to variations in the electric field near the surfaces of the detector, 
variations in the contaminants in the germanium, and not being 
able to sort out pure photoelectric events. 



B. Fan Source Scan 
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Fig. 6. A fan beam scanned across the side of the detector from the lithium 
side to the boron side. The x axis is the actual position of the source on the 
translation table and y axis is the depth of the interaction determined by taking 
the time difference in charge collection. The lower left hand corner corresponds 
to the front of the lithium face and the upper right hand corner to the front of the 
boron face. The bottom plot shows the location of the centroid for each position 
and a linear regression of the centroids. The error on the centroid is less than the 
diameter of the points. 



To test the depth resolution of the detector, the side of the 
detector was illuminated with a tightly collimated gamma-ray 
beam. A 1-mCi 57 Co source was mounted in a collimator con- 
sisting of two flat planes of tantalum approximately 1 1 .5 cm in 
length and 2 cm thick. The two planes are separated by 0. 1 -mm- 
thick spacers. This produces a well defined fan beam useful for 
scanning the detector. The fan source was scanned along the side 
of the detector using an x-y position table. The table has a posi- 
tion resolution of 0.025 mm and a range of 1 0.2 cm. The source 
was moved in 0.5-mm steps and data was collected at each point 
along the side of the detector. 

A histogram of the timing difference between charge collec- 
tion on each boron strip and any lithium strip was constructed. 
For each event, only one lithium and one boron strip could have 
a signal and their energies had to be within 5 keV of the 1 22 keV 



line. One boron strip in the middle of the detector was selected 
and the time difference for each position was plotted (see Fig. 6) . 

Based on a linear regression of the centroids for each position, 
the detector is shown to have an integral nonlinearity of 5.7% 
across the detector. This slight nonlinearity is probably due to 
changes in the electric field near the electrode structures on the 
faces. 

The time resolution of the system for the fan beam illumi- 
nating one position on the side of the detector is 14-ns full width 
half maximum (FWHM). This corresponds to 0.70 mm for this 
detector. The gamma-ray beam is 0.15 mm wide at the edge 
of the detector and the average electron motion at this energy 
is 0. 1 mm. Subtracting these contributions in quadrature from 
the overall resolution of 0.70 mm yields a depth resolution of 
0.68- ± 0.09-mm FWHM. 
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IV. Single-Detector Compton Telescope 

Having three-dimensional (3-D) readout of a germanium strip 
detector gives good position resolution in all three dimensions 
arid excellent energy resolution. This allows one to use a single 
detector as a Compton telescope, as opposed to the traditional 
configuration using two separate detectors in coincidence to 
measure two interactions. 

Consider gamma-rays coming from a point source. Some of 
these gamma-rays will Compton scatter in one location in the 
detector and then interact a second time at a second location in 
the same detector, depositing all of their energy in these two in- 
teractions. The Compton scattering angle in the first interaction 
can then be determined from the Compton formula 



cos 



where 0 is the Compton scattering angle, m e is the electron rest 
mass, and E\ and E 2 are the energies in keV deposited at the two 
interaction points. Knowing the position of the two interaction 
points can then be used to draw a cone of possible directions 
from which the gamma-ray source must be located. Drawing 
enough of these cones and determining the intersection point 
reconstructs an image of the gamma-ray source. 

This experiment was done with a 8.8 //.Ci 22 Na source placed 
41 cm from the boron side of the detector and a 1.3 ^Ci 137 Cs 
source 20 cm to the left of the 22 Na. The data was acquired 
for 45 minutes. Events that had two strips with signals on the 
boron side and two strips hit on the lithium side that added up 
to either 662 keV or 5 1 1 keV were selected. These events were 
then checked to make sure that each hit on the boron side had 
an exact energy match with a strip on the lithium side and that 
events were not in neighboring strips. This data set was then 
used to reconstruct the image using a simple ring sum algorithm. 

Each point on a plane located 41 cm from the detector was 
tested to see if it satisfied the Compton scattering formula within 
errors using the position and energy information from the event. 
Each pixel that satisfied these requirements was given a value 
weighted by the total number of pixels for each event. This was 
done for both orderings of the event since the true ordering is not 
always known. All events were then summed together, which 
produced the image shown in Fig. 7. A similar image was pro- 
duced when knowledge of the sources positions were used to 
determine the correct event ordering. 

Both sources are visible in the image and are separated by 
20 cm. The 137 Cs has better angular resolution because it has 
the higher gamma-ray energy. The position resolution is about 
5 cm, which corresponds to 7° angular resolution. This image 
would have been impossible without the depth resolution be- 
cause the interaction point would only have been defined by the 
overlapping front and back strips. 

V. Timing Methods 

All of the experiments in the previous section used the elec- 
tronics diagramed in Fig. 2. The TDCs were started by the CFD 
signals from our custom NIM boards and stopped by a delayed 
copy of the LED. To determine if the depth resolution mea- 
sured in the preceding sections is limited by the detector or the 
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Fig. 7. A reconstructed image of a 137 Cs and 22 Na source placed 41 cm from 
the boron face of the detector and separated by 20 cm. A Compton ring for each 
event was drawn at a distance of 41 cm from the front face of the detector and 
summed together to produce this image. 



electronics, the depth resolution was measured using different 
electronics setups. Commercial NIM modules were used to test 
these other timing methods due to their flexibility and ease of 
wiring. Due to channel limitations, only one lithium strip and 
three boron strips were instrumented. 

The CFDs implemented on our custom boards could be lim- 
iting the timing resolution of the system. To determine if this is 
the case, preamplifier signals were fed into Ortec timing filter 
amplifiers (TFA) set to 200 ns differentiation and integration 
time. The TFAs signal is sent to an Ortec CFD with 150 ns 
of external delay. The timing signals from the CFDs start the 
TDC channels, gate the ADCs, and, after being delayed, stop 
the TDC. To measure the timing resolution of this system, a 
57 Co fan beam illuminated a fixed position on the side of the 
detector which was 4.5 mm from the Uthium face. This beam 
illuminated all of the horizontal boron strips and, due to attenu- 
ation, the first few Uthium vertical strips. Using this electronics 
setup, the timing resolution was 9 ns, which corresponds to 
0.4 5 -mm depth resolution. Taking into account the beam width 
and electron motion, this system has a depth resolution of 0.41 
± 0.08 mm, which is better than the 0.68- db 0.09-mm resolution 
with the CFDs on our custom boards. The lower performance of 
the custom build electronics is probably due to jitter or noise in 
the design. Further tests will need to be done to determine the 
exact cause of the problem. 

Replacing the Ortec CFDs with LEDs and setting the trig- 
gering threshold to 20 keV resulted in a timing resolution of 9-ns 
FWHM as well. A comparison of the depth resolution for this 
configuration and for the configuration with the custom built 
NIM modules is shown in Fig. 8. 

The LED worked as well as the CFD for the 122 keV 
gamma -ray line but it is not known at this point if it would have 
the same resolution at a range of different energies. There are 
many other methods to determine timing accurately without 
the need for a delay line. One that we have implemented and 
will be testing soon uses a comparator to look at the crossing 
between the fast shaped preamplifier signal and the integral 
of the fast-shaped signal [7]. This timing circuit has been 
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F ig. 8. The solid curve is the time difference between charge collection on the 
boron and lithium sides of the detector using a LED to determine timing. The 
dashed curve is the time difference using the integrated NIM electronics used 
for the other experiments. 

produced in CMOS which would be useful for producing an 
ASIC that combines a shaped signal and timing information 
for an entire detector. All of these methods will be investigated 
further at a variety of energies. 

VI. Multiple Interactions 

As the gamma-ray energy increases, the likelihood of the 
gamma-ray depositing all of its energy in one pixel decreases. 
At a gamma-ray energy of 662 keV from a 137 Cs source, the 
photo-peak efficiency in the germanium detector is less than 
1%. The rest of the events will involve Compton scattering 
and charge sharing with neighboring strips. Depth information 
should be able to distinguish between these two types of events. 
This would allow charge sharing . events to be used in event 
reconstruction using the average position and the sum of the 
energies in the two strips. The Compton scattering events 
would then be available for reconstruction. This increases the 
efficiency of the detector by allowing more event types to be 
used in the final analysis. 

To look at depth information in neighboring strips, a L37 Cs 
source was placed 41 cm from the detectors boron side. Events 
were selected that had only one signal on the boron side and 
two neighboring strips hit on the lithium side. The time dif- 
ference between the two neighboring strips on the lithium side 
was histogramed. Charge sharing events should have essentially 
no time difference and Compton events should have a variety 
of time differences based on where the interactions occurred. 
Charge sharing events should be independent of source position 
while Compton scattering should be affected by source location 
because this movement causes changes in the scattering angles 
between the strips. This was all seen in the experiment as shown 
in Fig. 9. There is a center peak, charge sharing events, that 
was unaffected by source position and then Compton events that 
shifted with changing source position. More work is necessary 




-100 



-50 0 50 

Time Difference (nsec) 



Fig. 9. A 137 Cs source was used to illuminate the detector at two different 
points separated by 20 cm at a distance of 41 cm from the boron face of the 
detector. The time difference between neighboring strip hits on the lithium side 
is histogramed The solid line is for the source centered with the detector and 
the dashed line is for the source located 20 cm above the center. 



to make this technique useful for distinguishing between dif- 
ferent event types. 

VII. Conclusion 

The depth of a gamma-ray interaction can be measured in an 
orthogonal germanium strip detector to less than 0.5 mm. The 
depth information coupled with the x-A} position information 
from the strips yields a detector that is useful for a number of 
ground and space based instruments. Compton telescopes built 
from detectors with 3-D readout would have better image and 
energy reconstruction. Also, this enables the use of thicker de- 
tectors which would lead to less electronics for the same amount 
of detector volume. 
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Abstract 

We investigate two in situ calibration techniques (using 
gamma rays incident from the patient side) for PET detector 
modules that measure interaction depth via an analog ratio. 
For each crystal element in the camera, we need to determine 
two unknowns: the gain ratio K of the two detector (PMT and 
PD) signals and the depth dependence a of the signals (defined 
as the ratio of the signal observed when the crystal is excited 
at the end closest to and furthest from the detector). The light 
collection is assumed to be linearly dependent on distance from 
the detector end, in agreement with experimental results. The 
first method extracts K and tit using (a) the ratio of the PMT 
and PD signals for interactions that occur in the detector end 
closest to the patient (the most probable depth) and (b) the 
shape of the PMT pulse height distribution which reflects the 
exponential attenuation length in the detector. The second 
method utilizes the fact that H=PD+PMT (the total energy 
estimator) is independent of r=PD/(PD+PMT) (the depth 
estimator) when A' is correct, with a distribution position that 
is a dependent. On simulated data, both the gain ratio and 
depth dependence are determined with errors of 3% rms and 2% 
rms respectively, resulting in minimal degradation of energy 
and depth resolution. 

I. INTRODUCTION 

We have previously described the design for a PET detector 
module that consists of an array of LSO crystals coupled on 
one end to a photomultiplier tube and on the opposite end to 
an array of silicon photodiodes [1]. We calibrated this 
prototype detector module using an electronically collimated 
beam of annihilation photons incident from the side of the 
detector to excite at a specified depth of interaction (figure 1 a). 
Others have suggested calibrating different depth-encoding PET 
detectors using techniques similarly dependent on 
measurements with side incident beams to provide a photopeak 
pulse height versus depth look-up table [2-3]. However, using 
side incident gamma rays is awkward for complete cameras 
since recalibration is only possible by removing detectors 
from the gantry. 

We therefore investigate in situ calibration techniques {i.e. 
those that utilize only gamma rays incident from the patient 
side of the detector module such as in figure lb). We use a 
Monte Carlo simulation that includes the measured electronic 
noise, energy resolution, and signal levels, as well as energy 
dependent Compton and photoelectric cross sections in LSO. 



* This work was supported in part by the U.S. Department of 
Energy under Contract No. DE-AC03^76SF00098, in part by 
Public Health Service Grant Nos. P01-HL25840 and R01- 
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University of California Grant No. 1RB-0068. 



The light collection is assumed to be linearly dependent on the 
distance of interaction from the photodetector, which has been 
.empirically found to be a good approximation [I]. More 
precisely, the PD and PMT signals (in ADC counts) for an 
ideal detector are given as a function of depth x by: 

PDW = |v K«- 0 0 1) 0) 

PMT(*)=!* pmr [(cx- 11, (2) 

where k pc i and h pmt are the gain factors that convert energy 
deposition to ADC counts for each photodetector, E is the 
energy deposited in the crystal, L is the crystal length, and a 
is the depth dependence (defined as the ratio of the signal 
observed when the crystal is excited at the end closest to and 
furthest from the photodetector). Depth x is defined to be 0 at 
the PD end. Figure 2 shows these detector signals as a 
function of depth for an ideal detector, assuming noise free 
readout, perfect energy resolution, and monoenergetic energy 
deposition. 

The necessary calibration parameters are depth dependence 
(a) and the relative gain of the PMT and PD (K). More 
specifically, 



_ PMT(jc=L) _ PD(x^O) 
a ~PMT(x=0) PD(r=L). 



(3) 
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Figure I : Calibration geometry for a depth of interaction detector 
module. An electronically collimated beam of annihilation 
photons are incident from (a) the side of the detector module to 
excite at a specified depth of interaction or (b) the patient side of 
the detector module for in situ calibration. 
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Figure 2: The PD and PMT detector signals, which are linearly 
dependent on the depth of interaction (x) from the photodetector. 
The depth x is defined to be 0 at the PD photodetector. Total 
. energy is the sum of PD and PMT. 
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r _ _ PMTQr=L) PMT(*=ft) 
K ~ fa" PD(x=0) ~ a PD(x=0)' 



(4) 



In this paper we present two different //? situ calibration 
methods used to extract K and a for a single crystal element. 

n. Simulation 

Monte Carlo data are used to simulate the performance of a 
single LSO crystal coupled on one end to a silicon photodiode 
and on the opposite end to a photomultiplier tube. Normally 
incident 511 keV photons enter the crystal on the PD end and 
interact in the crystal via Compton scattering and/or 
photoelectric absorption, assuming energy dependent cross 
sections in LSO. A 1.2 cm attenuation length (l//i) of LSO is 
assumed [4]. Since there is only a single simulated crystal, 
Compton scattered photons that escape the crystal are lost and 
none are added from neighboring crystals. 

The photodetectors are assumed to have a quantum 
efficiency of 0.8 and 0.2 for the PD and PMT, respectively. 
The detected signal is smeared by a Gaussian distribution to 
simulate the effect of finite scintillator energy resolution. The 
rms width of the Gaussian is proportional to the square-root of 
the signal and the proportionality constant is chosen to match 
observed energy resolution in LSO. This scintillator based 
"noise" is summed in quadrature with the measured electronic 
noise ( 1 90 electrons rms for PD and 0 for PMT detector). The 
signal in terms of energy is multiplied by the gain factor k pc j 
(or k pmt ) to give the signal PD (or PMT) in terms of ADC 
bins. We want to extract the relative gain K and calibrate so 
that 1 PMT ADC bin equals I PD ADC bin after correction. 
With these assumptions, the total energy resolution at 511 
keV is -19% fwhm; this is pessimistic compared to the 
expected detector performance [1]. Calibration results extracted 
with a better energy resolution are found not to affect the 
results. 

in. Exponential Attenuation Method 

The light collection is linearly dependent on the depth. 
Thus, the exponential nature of annihilation photon 
interaction in the crystal is reflected as an exponential shape of 
the PMT and PD pulse height spectra. The first calibration 
method that we describe utilizes this exponential shape. 

Figure 3 demonstrates the underlying data that is used for 
the calibration. Consider 51 1 keV gammas that interact only 
via photoelectric effect in a noise free detector. Interactions at a 
single depth would appear as a delta function in the pulse 
height distribution. As the distribution of interaction depths is 
exponential (due to exponential attenuation in LSO), the PD 
and PMT pulse height spectra also have exponential slopes. 
The PMT spectrum has a clear peak at low pulse height bins, 
corresponding to interactions at the patient end of the 
scintillation crystal (x=0, the most probable depth). Because of 
light sharing, interactions at x~L {i.e. close to the PMT) 
populate high pulse height bins. Ideally, we could use the four 
edges of the spectra — PD(*=0), PD(x=Z,), PMT(x=0) and 



PMT(x=L) — to extract the calibration constants using 
equations (3) and (4). 

For a realistic detector that includes electronic noise, finite 
energy resolution, and CompLon scatter (as in figure 4), 'we 
can still identify the peaks corresponding to interactions at 
x=0. However, it becomes difficult to accurately identify the 
pulse height edges for the x=L interactions. Thus, we have 
developed an algorithm to determine K and a using the x-0 
peaks and the shape of the PMT spectrum. 

For a realistic detector, Gaussian smearing of the PMT 
pulse height spectrum rounds off the abrupt edges seen at x=0 
and x=L in figure 3 and alters the exponential shape near these. 
x values. However, the shape is unaffected for 0. \L£x < 
0.9L Thus the PMT pulse height spectrum is tit with an 
exponential over a limited region — from 5% above the x=0 
peak to 25% below the bin where the PMT spectrum tail first 
goes to zero — in order to extract the slope m (N =Nq txp(-b- 
mPMT) = Nq exp(-x/fl)). The ?MT(x=0) peak is estimated by 
the ADC bin with the maximum number of events. The 
combination of the peak PMT(x=0) and exponential slope m 
are used to extract alpha: 

L 

m fi PMT(x=0) ' 

The PD pulse height spectrum peaks for interactions 
corresponding to x=0, which is also estimated by the ADC bin 
with the maximum number of events. We can then obtain the 
relative gain factor K from PMT(x=fl) and PD(jc=0) using 
equation (4). 



a= 1 + • 



(5) 
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Figure 3: The PMT and PD distributions simulated with 0t=3 and 
K=\ for an ideal detector. Exponential attenuation (l//z) is evident 
for both PD and . PMT distributions, as well as the edges 
corresponding to x=0 and x=L. 
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Figure 4: The PMT and PD distributions simulated with a=3 and 
K=] for a realistic (solid lines) and ideal (dashed lines) detector. 
Exponential attenuation (l//i) is evident for the realistic PMT 
distribution (to the right of its peak), as well as the edges 
corresponding to x=0 . 
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We need to correct for systematic errors* such as those 
caused by the PMT(x=0) and PD(x=0) peak finding algorithm 
(see Figure 4). To obtain these correction factors, a data set of 
20 runs is simulated with different {a,K) combinations and 
200k events per run. The true (i.e. input) a and K constants 
range in value from 2-4 and 0.5-1.5 respectively to fully 
cover the expected detector conditions. The nominal values are 
ot=3 and K= 1 . This calibration data set is used to measure "raw 
values" of a and K with the exponential attenuation method 
described above. We plot the measured raw versus true values 
and fit a line for both ct and K, providing correction factors for 
future data sets. 

To estimate the errors in this method, a second data set of 
20 runs is simulated with 200k events each and the same 
(a,/Q combinations as in the calibration data set but with 
different random seeds. The exponential attenuation method is 
applied to each run, then corrected with the corresponding a or 
K calibration line. Figure 5 shows the resulting measurements 
for a and K as a function of the true values. A line of slope 1 
is also shown to represent the desired values. Thus, we obtain 
a with a 3% rms error and K with a 2% rms error using this 
method. In figure 6, the pulse height spectra and depth of 
interaction distributions are shown to be almost 
indistinguishable when using either the true or measured 
calibration constants. When exciting at a. fixed depth, we also 
see little difference in the depth distributions for true and 
measured calibration constants. Therefore, we have found that 
the calibration errors cause only a minute degradation of the 
energy and depth resolution. 

This technique uses only a crude estimation of the pulse 
height peaks (PMT(x=0) and PD(jc=0)) by using the bin with 
the maximum number of events for each photodetector. 
Clearly a more sophisticated fit for the peak region would give 
a better estimation. In addition, the PMT spectrum should be 
fit with an exponential after the Compton background has 
been subtracted. However, our investigation showed accurate 
slope extraction when fitting without subtraction. For speed 
and simplicity, we choose not to fit the peaks or subtract off 
the background since they aren't necessary. A run of 200k 
events is a very large sample when considering that we have 
over 20k crystals in need of calibration in our proposed PET 
camera. However, preliminary measurements (with 20k events 
per run) indicate that this technique will yield acceptable 
results even with a smaller number of events. Acquisition 
time for sufficient data* to perform calibration is approximately 
10 minutes (assuming 30k crystals, 20k events/crystal, and 1 
MHz acquisition rate). Computer analysis (histograming, 
fitting, and extracting constants) is expected to take 
approximately 1.5 hours on a typical UNIX workstation. 

IV. Constant Energy Sum Method 

Assuming the light collection to be linearly dependent on 
depth, the second calibration method utilizes the fact that the 
total energy estimator (E=PD+PMT) is independent of the 
depth estimator (f=PD/(PD+PMT)) if the relative gain K is 



correct. Figure 2 shows the PD, PMT, and total energy 
detector signals as a function of true interaction depth for an 
ideal detector with a correct relative gain of 1. Figure 7 shows 
the total energy as a function of the depth estimator for an 
ideal detector with either K=l (figure 7a) or K=\ .5 (figure 7b). 
A clear dependence on T is observed when the gain is 
incorrect, whereas the total energy is flat when the gain is 
correct. Thus wc can estimate K from the slope of E as a 
function of T. Figure 7 also shows that the total energy 
position is a function of the calibration constants; the line is 
centered for correct gain (figure 7a) and shifts to the left for 
K>\ (figure 7b). More precisely, the upper edge is given by 

n~o)=^. ( 6) 

Thus we can extract the calibration constant a by measuring 
the upper edge T(x=0) once we have determined K. 
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Figure 5: The correlation between true and measured calibration 
constants using the exponential attenuation method for (a) a and 
(b) K. The lines of slope 1 represent the desired values. 
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Figure 6: The (a) total energy pulse height spectrum and (b) depth 
of interaction distribution for both the true and measured 
calibration constants. The curves are almost indistinguishable. 
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Figure 7: The total energy estimator E=PD+PMT as a function of 
the depth estimator f=PD/(PD+PMT) for an ideal detector with 
ot=3 and (a) a correct relative gain of K=l or (b) an incorrect gain 
of K-\.5. Total energy is independent of f if the gain is correct. 
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Figure 8: The total energy estimator E=PD+PMT as a function of 
the depth estimator r=PD/(PD+PMT) for a realistic detector with 
a=3 and fa) a correct relative gain of K-\ or (b) an incorrect gain 
of /C=l.5. Total energy is independent of T if the gain is correct. 
The crosses indicate photoelectric events and the dots are 
Compton scatter events, as determined by an expectation 
maximization technique. The white lines are a fit of only 
photoelectric events (used to extract K). 

For a realistic detector that includes electronic noise, finite 
energy resolution and Compton scatter, the lines in figure 7 
become a scatter plot of E verses T but still demonstrate 
independence of T if the gain is correct (figure 8). Figure 8 
also shows an excellent separation between photopeak and 
Compton scatter events for both correct and incorrect gains. In 
order to properly extract K from the scatter plot for a realistic 
detector, we select only the photopeak events using an 
expectation maximization technique [5]. We model the 
probability densities of the events in the variables r and E as a 
mixture of two Gaussian functions and classify each event 
according to Bayesian inferencing [6]. The white line in figure 
8b, corresponding to the average E as a function of r for 
photopeak events only, is a hyperbola. Thus rather than fitting 
E, we make a linear fit of the inverse of E. We extract K from 
the slope and intercept of this line, giving an excellent fit for 
photopeak events as seen in figure 8. We are able to observe a 



clear upper edge r(x=0) for all K values. This edge is 
determined by histograming T and estimating the peak by the 
bin with the maximum number of events. We can then 
determine a from both K and r(jc=0) using equation (6). 

The same calibration data set used to measure a and K 
with the exponential attenuation method (described previously) 
is used for the constant energy sum method. However, only 
the first 20k events per run are used because of the 
computation time required by the expectation maximization 
algorithm. There are systematic errors such as the position of 
the r(*=0) peak, but the calibration results show a correlation 
between the measured raw values and the true values. We plot 
the measured raw versus true values and fit a line for both a 
and K, providing correction factors for future data sets. 

In order to estimate the errors, the constant energy sum 
method is applied to the same test data set used for the 
exponential attenuation method but using only the first 20k 
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Figure 9: The correlation between true and measured calibration 
constants using the constant energy sum method for (a) a and (b) 
K. The lines of slope I represent the desired values. 
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Figure 10: The (a) total energy pulse height spectrum and (b) depth 
of interaction distribution for both the true and measured 
calibration constants. Excellent agreement is seen. 
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events. These measured raw values of a and K are then 
corrected with the corresponding a or K calibration line. 
Figure 9 shows the resulting measurements for a and K as a 
function of the true values. Thus, we obtain a with a 14% 
rms error and K with a 4% rms error using this method. 
Although the rms error for a is much larger for this method 
when compared to the previous method (with 2% rms for a), 
the calibration errors cause only a minimal degradation of the 
energy and depth resolution (figure 10). For the results 
reported here, automated Matlab code was used to perform the 
calculations and it took approximately 2 minutes to calibrate 
each crystal (i.e. 20k events). Since running loops within 
Matlab is known to be very slow, an automated version in C 
is expected to be significantly faster. We anticipate that the 
constant energy sum method will be ~2 times slower than the 
previous method (i.e. ~3 hours for 30k crystals). 

V. DISCUSSION 

These results are based on a single simulated crystal. One 
cannot truly validate a detector calibration method or select 
among competing options without having constructed a 
number of functioning detector modules, • since presently 
unmodeled effects (such as Compton scatter from adjacent 
crystals, deviations from linear dependence of the light 
collection on depth, electrical and optical crosstalk between 
detector elements, etc.) may need to be accounted for. 
However, the concepts presented herein appear promising and 
practical based on this relatively simple model. 

Although we have assumed the light collection to be 
linearly dependent on depth for these methods (in agreement 
with experiment), they could easily be modified for a different 
depth dependence as long as the functional form is known; this 
can be determined using measurements with side incident 
gamma rays. Calibration techniques are by necessity detector 
specific and so it is difficult to accurately assess whether these* 
methods can be used with other detector designs. However, 
there are several current designs [2,3,7,8] that measure depth 
with an analog signal (or signals) whose amplitude depends on 
depth. Requiring the depth distribution to be consistent with 
the known exponential distribution should be sufficient to 
calibrate those with one analog signal, while requiring an 
additional constraint (i.e. requiring the sum to be depth 
independent or identifying interactions that occur at the crystal 
boundaries) should be sufficient for those that employ a ratio. 



vi. Conclusions 

We investigate two in situ calibration techniques for PET 
detector modules that measure depth via an analog ratio. The 
first method uses the ratio of the detector signals at the most 
probable depth as well as the exponential shape of the PMT 
pulse height spectrum in order to extract the calibration 
constants. We obtain the depth dependence (a) with a 3% rms 
error and the relative detector gain (K) with a 2% rms error 
using this exponential attenuation method. In the second 
method we select photopeak events and fit the total energy 
versus depth estimator (f). shape to extract K\ we then 



determine a from both K and the upper edge T(x=0). We 
obtain the depth dependence with a 14% rms error and the 
relative detector gain with a 4% rms error using this constant 
energy sum method. Estimated calibration times are 2-3 hours 
for a 30,000 crystal tomograph, dominated by fitting time. 
The larger errors in the constant energy sum method may be 
due to the smaller number of events (20k vs. 200k). However, 
both calibration techniques cause only minimal degradation of 
the energy and depth resolution so they are both viable options 
for in situ calibration. 
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Abstract 

Mercuric iodide, a high-Z, wide band-gap semiconductor material, is a room temperature 
gamma-ray detector. Poor hole transport properties limit the detector performance. By 
parallel pulse processing schemes we are exploring this problem and making corrections to 
enhance resolution. This yields interaction depth information which is used to estimate hole 
collection deficits. Deficit estimates are used as corrections to pulse height to enhance 
resolution, or to reject pulses from high-deficit regions. These techniques have been 
implemented by several schemes for filtering the pulses using off-the-shelf amplifiers and 
computer simulation of amplifiers on detectors from 1-3 mm thick. Hardware and simulated 
results are compared and improved spectral performance is discussed. A resolution better 
than 7% for correction and 2.5% for rejection of hole deficit is obtained. 



*1 This work was performed under the auspices of the US Department of Energy under 
Contract No. DE-AC08-83NV 10282. 
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Spectroscopic-Grade X-Ray Imaging up to 100-kHz 
Frame Rate With Controlled-Drift Detectors 

Andrea Castoldi, Member, IEEE, Chiara Guazzoni, Member, IEEE, Pavel Rehak, Member, IEEE, and Lothar Striider 



Abstract— Controlled-drift detectors are fully depleted silicon 
detectors for X-ray imaging that combine good position resolu- 
tion with very fast frame readout The basic feature of the con- 
trolled-drift detector is the transport of the charge packets stored 
in each pixel column to the output electrode by means of a uni- 
form drift Held. The drift time of the charge packet identifies the 
pixel of incidence. Images of an X-ray source obtained with the 
con trolled-drift detector up to 100-kHz frame rate are presented 
and discussed. The achievable energy resolution as a function of 
the operating temperature and frame rate is analyzed. 

Index Terms — Controlled-drift detector, fast readout, X-ray 
imaging. 

I. Introduction 

THE controlled-drift detector (CDD) was proposed in 1997 
[ 1 ] , and the first experimental evidence of its working prin- 
ciple was reported in [2]. The device is built on a fully de- 
pleted 300-/xm-thick high-resistivity wafer and is operated in in- 
tegrate-readout mode. The basic idea of the CDD is to generate 
columns of equally spaced potential wells for the electrons by 
superposing a periodic perturbation of sufficient amplitude to 
a linear drift potential. During the integration mode, the signal 
electrons are stored within these wells. The removal of these 
potential barriers in an externally controlled way allows the use 
of the linear drift potential to sweep the electrons toward the 
readout electrodes. The time between the removal of the barriers 
and the arrival of the signal electrons at the readout electrodes 
gives the position of the irradiated pixel along the drift direc- 
tion. The second coordinate is obtained by providing a separate 
readout electrode for each pixel column. Arrays of deep /^im- 
plants provide lateral confinement. In Fig. 1 , the layout of the 
anode region of the controlled-drift detector is shown. 

This readout mechanism is inherently faster than the clocked 
transfer of the rows of pixels toward the readout section typical 
of the pn charge-coupled device (CCD) because now charge 
transfer and signal processing are carried out simultaneously. 
The frame rate in the CDD is ultimately limited by the 
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Fig. 1 . Layout of the anode region of the controlled-drift detector. The pixel 
size is also indicated. 



electron drift time, which is typically 2-6 /us for 1 -cm-long 
device depending on the applied drift field. The advantage 
of this readout mechanism is twofold: 1) a higher frame 
rate with respect to the fully depleted pn-CCD [3], which 
represents the reference device for high-resolution imaging 
and spectroscopy of X rays, and 2) a lower contribution of 
the thermal noise due to a shorter integration time, leading 
to high energy resolution. 

We performed the first extended test of the imaging and spec- 
troscopic capability of the CDD by irradiating a column of the 
detector at room temperature. The aim of the measurements was 
twofold: 1 ) verify the proper operation of the CDD at frame fre- 
quencies close to the limiting value set by the electron drift time 
and 2) verify the expected improvement of the energy resolution 
with the frame frequency. In Section II, we present the images 
of a 55 Fe source obtained with the CDD at different frame rates 
up to 100 kHz. In Section III, the achievable energy resolution 
is discussed as a function of the frame rate and of the operating 
temperature. 

II. X-RAY IMAGING 

One column of the detector has been irradiated with a 00 Fe 
source to verify the imaging capability, the maximum frame 
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rate, and the achievable energy resolution. The pulses delivered 
by the on-chip JFET (source follower configuration) were fed to 
a low-noise voltage preamplifier followed by a pseudo-Gaussian 
shaper (0.25 -/xs shaping time) and an 8-bit digitizer controlled 
by a PC to acquire the output waveforms, subtract the back- 
ground (mainly due to thermal generation), and compute the 
amplitude and the drift time of each pulse. We acquired several 
images of a 55 Fe source at room temperature for different frame 
rates. As the detector area is unshielded, the ratio between the 
integration time and the readout time is typically chosen greater 
than one in order to limit the number of out-of-time events. In 
our case, it was set equal to nine in all measurements. With this 
choice, we tested the CDD up to a 100-kHz frame rate, corre- 
sponding to a duration of the readout phase of 1 lis. Considering 
that the electron drift occurs within the first 0.7 /las of the readout 
phase, for the 100-kHz case signal processing was overlapping 
with the subsequent integration phase. 

The operating drift field was 200 V/cm and the amplitude of 
the surface perturbation was 2 V. 

Fig. 2(a) shows the scatter plot energy versus drift time (i.e., 
position) of the detected events when the CDD is operated at a 
1 0-kHz frame rate. Moving along the time axis, we see that the 
events are gathered in well-separated clusters centered at about 
6 keV (the Mn Ka and Mn K/3 lines are not resolved), indicating 
the illuminated pixels. In Fig. 2(b), the scatter plot related to a 
30-kHz frame rate operation is reported. As we can see, the clus- 
ters relative to the K/5 line start to separate from the ones corre- 
sponding to the Ka line as the energy resolution is improving. 
In the image taken at a 100-kHz frame rate [Fig. 2(c)], the Ka 
and K/i lines can be clearly distinguished. 

A more detailed analysis is obtained by plotting the measured 
events against one parameter at a time. Fig. 3 shows the distribu- 
tion of the events along the time axis for the image acquired at a 
1 00-kHz frame rate. The full-width at half-maximum (FWHM) 
of the peaks is smaller than 1 5%, leaving a margin for reduction 
of the pixel size along the drift coordinate (actual size is 180 
/mi). 

Fig. 4 shows the total distribution of the event energies (i.e., 
the spectrum of the 55 Fe source collected by all the pixels) for 
two different frame rates (10 and 100 kHz). 

At a 10-kHz frame rate, the energy resolution at the Mn Ka 
line (5.899 keV) is 837 eV FWHM, corresponding to an equiv- 
alent noise charge (ENC) of 98 electron rms. By increasing the 
frame frequency, we shorten the integration time and therefore 
reduce the amount of leakage charge that accumulates in the 
pixels, thus improving the energy resolution. The energy reso- 
lution at 100-kHz frame rate is 339 eV FWHM, corresponding 
to an ENC of 37 electrons rms. 

To further reduce the integration time, we changed the ratio 
between the duration of the integration and the readout phase. 
In Fig. 5, the scatter plot energy versus drift time of the detected 
events when the CDD is operated at a 1 00-kHz frame rate, with 
a duration of the integration phase of only 6.7 p-s. As expected, 
we obtained a better energy resolution at the Mn Ka line equal 
to 296 eV FWHM, corresponding to an ENC of 32 electron rms 
(see Fig. 6). 
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Fig. 2. Scatter plot energy versus drift time of the X-rays of a "Fe source 
collected by a column of the CDD. The measurement was carried out at room 
temperature. The clusters have been labeled with the corresponding pixel 
numbers, (a) 1 0-kHz frame rate (integration time 90 ^s), (b) 30-kHz frame rate 
(integration time 30 /is), and (c) 100-kHz frame rate (integration time 9 /<s). 
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III. ANALYSIS OF THE ENERGY RESOLUTION 
The readout noise can be expressed in terms of the ENC as 

ENC 2 = ENC 2 ! + ENC 2 ^ + ENC^ ig . (1) 

The noise analysis shows that the contributions are from the 
noise of the electronic chain (ENC cl = 20 electrons rms), the 
quantization noise of the 8-bit digitizer (ENC dig = 12 elec- 
trons rms), and from the parallel noise (ENCparaifei) due to the 
leakage current (about 40 pA/channel) integrated in the pixels 
during the integration time. 

We can model the noise contribution due to the integrated 
leakage charge as follows. During the integration mode, the 
leakage current fills the potential wells of the pixels. The av- 
erage value of the leakage charge in each pixel is equal to 





8000 




7000 












6000 


s 










5000 




4000 




3000 



0.0 



0.1 0.2 0.3 0.4 
time [^s] 



0.5 



0.6 



Fig 5 Scatter plot energy versus drift time of the X-rays of a 55 Fe source 
collected by a column of the CDD at 1 OWcHz frame frequency. The integration 
time was set to 6,7 /ts and the readout time to 3.3 ps. 
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where 

J L leakage current per unit area; 

/lpixel area of one pixel; 

r int duration of the integration time. 
After the transition to the drift mode, the charges Q L are 
released by the pixels and start moving toward the collecting 
anode. The motion of the integrated charge pulses gives rise to 
a current pulse that lasts for a time equal to the drift time from 
the last pixel (T driftma J and of average amplitude 

J,iUi 2drift milJl Adrift 
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2drift„ 



(3) 



Ql, — Jl ' -^pixel ' ^int 
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where JV pix «i is the number of pixels in a column and 1 L is 
the leakage current at the anode when the CDD is operated in 
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Fig. 7. Measured ENC 2 as a function of the inverse of the frame rate. The 
derived contribution of the parallel noise is also shown. 

free-running. The last equality is a slight overestimation due to 
the fact that 1 L also includes the leakage current generated in the 
anode region. As the drift time from the last pixel (7drift mftX — 
0.7 lis) is greater than the shaping time (0.25 ^s), we can esti- 
mate the parallel noise contribution (ENC£ aPalleI ) to the ENC 2 
as follows: 



ENC 2 )arallel - A 3 qI L 



1 + 



(4) 



where 

A3 shape factor of the filter for parallel noise contribution; 

q electron charge; 

r s }j time constant of the filter. 
A first confirmation of the proposed model is given by com- 
paring the energy resolution at the Mn Ka: line of the events . 
collected during the integration phase with one of the so-called 
out-of-time events (i.e., the events hitting the detector during the 
readout phase and whose position of incidence cannot be de- 
termined). The two spectra are compared in Fig. 6 for a frame 
frequency of 100 kHz and a duration of the integration phase 
of 6.7 /us. We considered only the out-of-time events hitting the 
detector in the last 0.3 /xs of the readout phase when the leakage 
current pulse of amplitude Iunt has expired. The resolution at 
the Ka line is 237 eV FWHM equal to the resolution achievable 
with the controlled-drift detector operated in free-running, that 
is, always in the readout phase. 

In Fig. 7, the measured ENC 2 is reported as a function of the 
inverse of the frame rate. Subtracting quadratically the known 
contributions of the noise of the electronic chain and the quanti- 
zation noise from the measured ENC, we calculated the parallel 
noise contribution. At 1 0-kHz frame rate, the parallel noise con- 
tribution is of 95 electrons, while at 1 00-kHz frame rate the par- 
allel noise contribution reduces to 3 1 electrons. From the slope 
of the fit, we get II = 40 pA. The plot of ENC 2 as a function of 
the frame frequency shows good agreement between model and 
experiment. The spread of the measured points is also related to 
the approximated energy-voltage calibration factor. 

As verified, the increase in the frame rate allows us to reduce 
significantly the dominant noise contribution of the integrated 
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Fig. 8. Contour plot of the achievable energy resolution in terms of ENC as a 
function of the detector working temperature and of the frame frequency. The 
ratio between integration time and readout time is set equal to nine. 

leakage charge. A well-known way to reduce the thermally gen- 
erated current is to cool down the detector. In the case of the 
fully depleted pn-CCDs [3], the only way to reach maximum 
performances is to cool down the detector, as the duration of the 
readout time is typically limited to the millisecond range. 

The achievable energy resolution, expressed by (1), has been 
estimated as a function of the operating temperature and of 
the frame rate. To this end, we assumed that the temperature 
dependence of the detector leakage current is the same as n?(T) 
(where ni(T) is the intrinsic carrier density at temperature 'I 7 ) 
while the noise contribution of the front-end electronics and of 
the digitizer was considered constant. Neglecting the increase 
of the JFET trasconductance at lower temperatures leads to a 
pessimistic estimate of the energy resolution. The ratio between 
integration time and readout time is set equal to nine. The 
equire solution lines of the estimated ENC in the plane (frame 
frequency, temperature) are shown in Fig. 8. 

IV. Conclusion 

In this paper, we presented the first extended characteriza- 
tion of the position-sensing and spectroscopic capabilities of the 
CDD. One-dimensional images of a 55 Fe have been acquired 
with the CDD operated at room temperature at different frame 
frequencies up to 1 00 kHz. The achieved energy resolution at the 
Mn Ka line and at 1 00-kHz frame frequency is 339 e V FWHM. 
The detailed noise analysis that has been presented shows that 
a dominant contribution to the energy resolution is that of par- 
allel noise. Moreover, the use of a digitizer with a number of bits 
greater than eight will practically cancel the quantization noise 
contribution. The possibility to operate the detector at higher 
frame rates (i.e., shorter integration times) is equivalent to virtu- 
ally cooling the detector. As a result, the room-temperature en- 
ergy resolution approaches the limiting value, obtainable only 
at cryogenic temperatures if the detector is operated at lower 
frame rates. 

We can expect to achieve energy resolutions close to 
present state-of-the-art X-ray imagers with more compact 
Peltier-cooled controlled-drift detector systems. 
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Abstract 

Hgl 2 detectors fabricated with pixelated anodes were tested as potential room temperature spectrometers. Spectra 
were obtained with 5 and 10 mm thick devices, and their results were compared. The 5 mm thick devices showed some 
effect of the anode weighting potential, while the 10 mm thick devices showed the effects of charge diffusion/sharing and 
electron trapping. A resolution near 1.4% FWHM using a 137 Cs source was achieved after correcting for the interaction 
depth within the detector. The resolution from individual depths was as low as approximately 1.1% FWHM. 
© 2003 Elsevier Science B.V. All rights reserved. 
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1. Introduction 

Mercuric iodide (Hgl 2 ) crystals are being 
investigated as a potential y-ray spectrometer that 
can operate at room temperature. However, like 
most compound semiconductor y-ray spectro- 
meter materials, Hgl 2 suffers from poor charge 
transport characteristics. Despite the fact that 
Hgl 2 has a wide band gap, high average atomic 
number, and a high density, the crystals suffer 
from low electron and hole mobility, significant 
hole trapping, and material non-uniformity. For 
these reasons, Hgl 2 spectrometers have been 
limited to a thickness less than 3 mm. To improve 
the spectroscopic performance of thicker Hgl 2 
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detectors, one must reduce the depth dependence 
of the generated pulse from a y-ray interaction 
event common to these materials. 

The use of single-polarity charge sensing was 
implemented in order to overcome the severe hole 
trapping present in wide band gap semiconductor 
detectors, including Hgl 2 . The introduction of 
pixelated anodes to CdZnTe [1,2] and Hgl 2 [3] 
showed that the resolution of wide band gap 
semiconductor detectors can be dramatically im- 
proved compared to standard planar electrode 
configuration. By simultaneously measuring the 
signal from the planar cathode and anode pixels, 
the interaction depth [4] can be obtained by 
calculating the ratio of the induced charge on 
each electrode. 

In this paper, we present and compare spectro- 
scopic results from 5 and 10 mm thick Hgl 2 
detectors. Both raw and depth-corrected spectra 

rights reserved. 



192 



J.E. Baciak, Z. He I Nuclear Instruments and Methods in Physics Research A 505 (2003) 191-194 



are presented, along with a comparison of the peak 
position as a function of depth of interaction. 



2. Detector description and setup 

The anode pixel configuration used in this work 
is shown in Fig. 1. The four small pixels are 
surrounded on all sides by a large anode plate. The 
area of each pixel is ~ 1 mm 2 and is comprised of 
palladium. When the cathode is negatively biased, 
electrons moving toward the anode surface induce 
charge on the pixels. The induced charge on a pixel 
is determined mainly by the movement and 
collection of electrons underneath each pixel. The 
signal produced by the anode pixel is only slightly 
dependent upon the distance traveled by the 
electrons, and is not significantly affected by hole 
movement. 
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Fig. 1. Top: view of the Hgl 2 anode configuration. Bottom: 
cross-sectional view of the crystal. The lines under the pixels 
represent the region of the detector where the majority of the 
induced charge on the anode is produced. 



The detector electrodes are connected to Amp- 
Tek [5] A250 charge-sensitive preamplifiers, which 
are connected to standard N1M shaping amplifiers 
for pulse shaping and amplification purposes. The 
signals are fed into peak-hold circuitry to hold the 
peak of the shaped pulse for a sufficiently long 
time before the ADC samples the pulse amplitude. 
The signals are analyzed by programs produced 
using LabVIEW™ and Matlab® software. 



3. Results 

3.1. 5 mm thick detectors 

The majority of the data from 5 to 6mm thick 
detectors have been described previously [6]. These 
detectors had resolutions near 2.1% FWHM at 
662 keV after correcting for the interaction depth. 
Interaction depth is broken up into 20 thin strips 
known as the "depth index." For example, a depth 
index of 1 would indicate events near the anode 
surface, while a depth index of 20 indicates events 
near the cathode surface. The resolution obtained 
at individual depth indices reached 1.6% FWHM 
(Fig. 2). The resolution of the spectra was inde- 
pendent of the depth index; however, spectral 
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Fig. 2. 137 Cs spectrum collected over 22 h from an individual 
depth of a 5.9 mm thick detector (Pixel #2) biased at -1500 V 
and using a shaping time of 4 us. The resolution of the 662 keV 
y-ray was 1.55%. The events selected here are from interactions 
near the middle of the detector. 
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resolution from events near the electrodes (anode 
or cathode) degraded to 2.3-3.7% FWHM. Poor 
resolution near the anode was expected because, in 
this region, the induced charge (weighting poten- 
tial) has a steep slope. Thus, the recorded pulse can 
be highly dependent on the interaction depth near 
the anode and on the induced charge produced by 
the small amount of hole movement present in 
these detectors. For events near the cathode, the 
degradation in energy resolution was due to 
electron trapping, charge sharing between the 
anode pixel and the non-collecting anode, and 
material non-uniformity effects. 



3.2. 10 mm thick detectors 

Fig. 3 is an example spectrum collected over the 
same time frame as the 5 mm devices. The 
resolution of the spectrum was approximately 
3% FWHM at 662 keV prior to correcting for 
interaction depth, and improved to 1.4% FWHM 
after depth correction. Resolutions as low as 1.0- 
1.1% FWHM can be observed when analyzing the 
spectra from individual depths (Fig. 4). The lack of 
dependence on hole movement and the depen- 
dence of peak position on interaction depth were 
also seen in the 10 mm detector. 
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Fig. 3. 137 Cs spectrum collected over 22 h from a 10 mm thick 
detector (Pixel #2) biased at -2500 V and using a shaping time 
of 8 us. The electronic noise of the system was approximately 
2.8 keV. This value was consistent for all pixels tested. 
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Fig. 4. l37 Cs spectrum collected over 22 h from an individual 
depth of a 10 mm thick detector (Pixel #2) biased at -2500 V 
and using a shaping time of 4 us. The resolution of the 662 keV 
7-ray was 1.06%. The depth index of this spectra was 10, which 
corresponds to v-ray interactions from the middle of the 
detector (i.e., 5mm from the cathode or anode surface). 



4. Discussion 

It is important to note that differences in the 
weighting potential, electron trapping, and charge 
sharing between anode pixels can have a signifi- 
cant effect on the spectra between 5 and 10 mm 
thick Hgl 2 detectors. For thinner pixelated detec- 
tors, there will be a larger dependence on depth as 
electrons traverse the detector. This is shown in 
Fig. 5. The slope of the weighting potential 
throughout the majority of the thickness of the 
detector is slightly larger for 5 mm detectors. The 
non-zero slope in the bulk of the detector causes 
the induced signal to have a slight depth depen- 
dence This results in a variation of photopeak 
position as a function of interaction depth. This 
effect can be seen when comparing peak position 
data with expected results (Fig. 6). 

When comparing the peak position as a function 
of depth index between 5 and 10 mm thick 
detectors, a number of features can be distin- 
guished. The 5 mm peak positions tended to lie 
below the expected curve for nearly all depth 
indices. This indicates that there was a more severe 
depth dependence on the induced charge than 
what was expected. The assumptions on weighting 
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Fig. 5. Weighting potential versus normalized depth for 5 and 
10 mm thick Hgl 2 pixelated anodes. The weighting potential for 
the planar cathode is also presented for reference. 
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Fig. 6. Normalized peak position as a function of depth index. 

potential (such as material uniformity) may not be 
entirely correct in this case [7]. A more severe 
depth dependence would adversely affect the 
overall resolution even with depth correction. 
Despite its advantages, depth correction cannot 
overcome the finite difference in the weighting 
potential within an individual depth (or voxel) 
itself. 

The number of depths that contained a clear 
photopeak was larger for 10 mm thick detector as 
a result of the weighting potential While the 10 mm 
measurements agree quite well with expectations, 
the measured data tends to be slightly higher for 
events near the anode (depth index 3-7) and 
slightly lower than expected for events near the 



cathode (depth index 15-20). The slightly higher 
induced charge could be due to a small amount of 
hole movement within the detector. The decrease 
in induced charge from y-ray events near the 
cathode could be due to electron trapping and 
charge sharing. 



5. Conclusion 

Hgl 2 detectors of 5-10 mm thickness with 
pixelated anodes have been tested as potential 
room temperature y-ray spectrometers. Since hole 
collection is unnecessary, relatively low electric 
fields and short shaping times were used. By using 
pixelated anodes and depth sensing and correction 
techniques, the resolution of individual pixels can 
be less than 2% FWHM when using a 137 Cs 
source. These detectors show no significant polar- 
ization effects resulting in the degradation of 
resolution. However, the 10 mm results showed 
the effects of charge trapping and charge sharing 
that degraded the signal from events near the 
cathode. If these effects can be reduced, the 
detector performance (resolution and photopeak 
efficiency) should improve. 
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